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Abstract 

We present results of a decay width calculation for a hybrid exotic meson(/i, J = 1 in 
the decay channel h — > Tva%. This calculation uses quenched lattice QCD and Liischer's finite box 
method. Operators for the h and 7ra± states are used in a correlation matrix which was expanded 
by varying the smearing and fuzzing levels at source and sink points. Scattering phase shifts for 
a discrete set of relative irai momenta are determined using eigenvalues of the correlation matrix 
and formulae derived by Liischer. The phase shift data is very sparse, but fits to Breit-Wigner 
models are made, resulting in a decay width of about 80 MeV. 
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INTRODUCTION 



Most states in the known hadron spectrum are resonances. Their decay widths, arguably 
more so than their masses, are sensitive indicators of the intrinsic structure and forces 
between the constituents. The decay of molecule-like states, for example, may be driven by 
a residual (effective) force and thus give rise to widths typical of nuclear physics. Spatially 
more compact hadrons could follow a quite different mechanism. 

Decays can be studied with lattice techniques [1]. In particular, lattice QCD provides 
us with ab initio access to unusual hadrons, like hybrids and exotics, which are a focus 
of some experimental programs currently under way. We here give an overview about a 
quenched lattice simulation set up to estimate the decay width of a J PC = 1 h hybrid 
exotic meson resonance, say h, coupled to a nai decay channel. Correlation matrices based 
on the corresponding operators yield (excited) mass spectra. Liischer's finite box method 

SB 

is then applied to get the scattering phase shifts at a discrete set of relative nai 
momenta. A decay width is determined by fitting Breit-Wigner functions. 

Such a simulation is viable, on a reasonably sized lattice, if the energy spectrum exhibits a 
level crossing between states excited by h and operators, preferably, close to the physical 
value of the pion mass. Unlike in p — > tttt, where this is very hard to achieve, h — > irai allows 
a relative s-wave in the decay channel and may thus proceed with relative lattice momentum 
zero. The feasibility of the current work much relies on this fact. 



LATTICE SIMULATION 



Hybrid mesons are quark-antiquark pairs having valence gluons as structural components. 



Following 



our choice for the hybrid operator is 



w = iEEiMwM^M - ftftft)). (i) 



1=1 



Here a, b are color indices and Fij is a product of SU(3) link matrices defined on paths 
forming a clover leaf in the i— j plane with center at x, and V is the spatial lattice volume. 
The quantum numbers are 1 = 1 and J PC = where the combination F^ — Fj- is needed 
for positive charge conjugation. Appropriate operators for nai can be constructed from 

= —= Yl 6 y - x,f d{xt)^u{xt) J(yt)7 5 7 i rf(yt) . (2) 

x y 



The simplest spatial configuration for an s-wave {A\ irrep) is the choice r = 0, which we use 
in this work. The elements in the resulting correlation matrix, 

C(*,tb)= ( ° h ' h{t,t0) Ch ^ to) \ , ( 3 ) 

are understood to include spin traces, for example Ch,h(t, t ) = Y^=i(Oh+;j(t)Of l+ .Ato)), and 
so on. 

The flavor structure of the operator (j2J) is ixai ~ dudd, which gives rise to (two) equal 



time contractions d(t)d(t) in Q. To avoid having to compute those propagator elements we 
observe that these do not occur for ttKx ~ duds. Since the experimental a\ and Ki masses 
are very close, 1260 MeV and 1270 MeV respectively, and given that only the mass spectrum 
is input to Liischer's method we argue that neglecting equal time contractions in C nai:7Tai 
should not significantly effect our results. It should be emphasized that this approximation 
only is applied to C 7raij7rai in (j3J), the off-diagonal correlators can be computed or inferred 
using hermiticity. 

Wuppertal style smearing is done on the quark fields, using APE style gauge link fuzzing 
in the process. We use a common strength factor (a = 2.5) and either 1, 2, and 3 smearing 
iterations. The same smearing is done at source and sink points. In this way the correlator 
(J3J) expands to a 6 x 6 hermitian matrix 

C(t,t )2x2 > C(t,t ) exe . (4) 

Simulations are made using the anisotropic Wilson gauge action with bare aspect ra- 
tio a s /a t = 2 on lattice sizes L 3 x 24, L = 8,10,12, at (3 = 6.15, and the corre- 
sponding Wilson fermion action in quenched approximation. Four hopping parameters, 
k = 0.140, 0.136, 0.132, 0.128, were used with a multiple mass solver [5J]. Figure. [T] shows 
the ground state masses of the h and the nai systems versus m 2 ^. These are obtained using 
effective mass functions based on output of individual h, ir, and ai correlators at smearing 
level two. The three parameter fit model, see inset in Fig. Q matches most of the energy- 
versus-m^ data produced in this work. Its choice is purely empirical. Setting the physical 
scale to the p meson mass results in a lattice constant of a t = 0.07 fm, and this value is used 
consistently for the width calculations. 
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FIG. 1: Masses of the h and tt + a\ non-interacting system versus the pion mass squared. The 
extrapolation shows an energy level crossing at light pion mass. If the physical scale is set to the 
ai meson mass, 1.26 GeV, then the hybrid mass extrapolates to 1.73 ± .39 GeV. If the p meson 
sets the scale, then the hybrid mass extrapolates to 1.46 ± .31 GeV. 

ANALYSIS 

To extract the mass spectrum the 6x6 correlation matrix (£Q) is first diagonalized on a 
fixed timeslice ti, 

C(ti, t ) = V(h, t )£>(ti, t )V^(t l7 1 ) , (5) 

where V(ti, to) is unitary and D(t\, to) is diagonal and positive definite. To ensure the latter 
t\ should be an 'early' timeslice, we here use t% — to = 4. The correlation matrix is then 
subject to a basis transformation and normalization according to 

1 .... . 1 



A(t) 



:y t (tl,to)C(t,t )\/(tl,t )- 



(6) 



/£>(Mo) ^(ti.to) 
This procedure amounts to solving the usual generalized eigenvalue problem. The only 
modification being that the eigenvectors (columns of V) used for the basis transformation in 
(JBJ) are employed at a fixed timeslice. Thus we find that A(t) is approximately diagonal. We 
will simply refer to the diagonal elements, A nn (t), as the eigenvalues, A n (t). Compared to 
'diagonalizing at all t' statistical fluctuations are reduced. The eigenvalues A n (t), n = 1 . . . 6, 
give rise to effective mass functions from which the spectral energies W n are obtained. Energy 
spectra for the 12 3 x 24 and 10 3 x 24 lattices are shown in Fig. 121 
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FIG. 2: Energy spectra W n , n = 1 . . . 6, and extrapolation to small pion mass for the 12 3 x 24 and 
10 3 x 24 lattices. Error bars are omitted for clarity. 

With the spectra known, subsets of masses W n are selected that fall within the elastic 
region, 

(m w + m ai ) < W n < 2(m n + m ai ) . (7) 

All masses in (JJJ) are taken at each of the four values of m v and in the limit m n — > using 
the extrapolation. For those sets the relativistic dispersion relation 



W n = y/ml + k 2 n + Jm 2 ai + kl 



(8) 



is solved for the relative nai momenta k n . Note, k n is a prediction for the relative momentum 
of a irai system and is not subject to lattice discretization. These momenta are input to 
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Liischer's formula [2| for the s-wave scattering phase shifts S n , 

. Tr 3 ^ 2 q n k n L s 
tano n = — — rr , q n = • (9) 

Here 2(1; q 2 ) is a generalized ^-function, and L s = La s is the physical size of the spatial 
box (using the bare anisotropy, a s = 2a t ). We then attempt fits to all sets of phase shift 
data with a Breit-Wigner model, 

tana(fc) = jfrj^FTM > with W ^ k ) = V m l + k " + + p > ( 10 ) 

where i?o and T are parameters. When the fit is successful, T may be interpreted as a 
decay width. In Fig. El we show selected results for the scattering phase shift data and the 
corresponding Breit-Wigner fits. It is apparent that most of the scattering phase angle data 
does not alone justify the use of Breit-Wigner functions. Use of these functions requires one 
to assume apriori that a resonant condition exists. 



RESULT 



Physical decay widths, found by using extrapolated energy spectra, varied from 78 ± 28 
MeV to 88 ± 23 MeV for the three lattice volumes tested. A large portion of the statistical 
error came from the hybrid correlator. As shown in Fig. El the phase shift data is very 
sparse. The best fits to the Breit-Wigner model returned a resonance energy of E « 1.9 
GeV. 

Phase shift data obtained using non-extrapolated spectra did not match up well to the 
Breit-Wigner model. The reason being, as shown in Fig. ^ the pion masses used are quite 
far from the actual energy level crossing, and thus, one is less likely to observe resonance 
behavior. 



SUMMARY 

Using Liischer's method, a decay width for the 1~ + exotic meson is calculated to be 
about 80 MeV with a statistical error of 25 MeV. The number of data points available to 
fit Breit-Wigner functions is sparse. There is a definite energy level crossing between the h 
and iTdi system at light pion mass giving credence to the simulation. 



6 



180 
150 
120 
90 
60 
30 


90 
75 
60 
45 
30 
15 




1 I t I 1 I 1 








/ r-. 088(23) GeV 




/ 0=6.15 




12 3 x24 






1 1 1 1 



0.0 



0.0 



0.2 



0.2 



0.4 0.6 
k [GeV] 



0.8 



1.0 




r = . 078(28) GeV 



(3=6.15 

10 3 x24 



0.4 

k [GeV] 



0.6 



0.8 



FIG. 3: Scattering phase shift versus momentum for the 12 3 x 24 and 10 3 x 24 lattices using 
extrapolated spectra from Fig. The curve fits are Breit-Wigner models (jlOj) . giving widths of 
88 ± 23 MeV and 78 ± 28 MeV respectively. 
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